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ABSTRACT 

The scattering polarization in the solar spectrum is traditionally modeled with each spectral line treated separately, but this is generally 
inadequate for multiplets where ./-state interference plays a significant role. Through simultaneous observations of all the 3 lines of a 
Cn triplet, combined with realistic radiative transfer modeling of the data, we show that it is necessary to include ./-state interference 
consistently when modeling lines with partially interacting fine structure components. Polarized line formation theory that includes 
./-state interference effects together with partial frequency redistribution for a two-term atom is used to model the observations. 
Collisional frequency redistribution is also accounted for. We show that the resonance polarization in the Cn triplet is strongly 
affected by the partial frequency redistribution effects in the line core and near wing peaks. The Cr i triplet is quite sensitive to the 
temperature structure of the photospheric layers. Our complete frequency redistribution calculations in semi-empirical models of the 
solar atmosphere cannot reproduce the observed near wing polarization or the cross-over of the Stokes Q/I line polarization about 
the continuum polarization level that is due to the ./-state interference. When however partial frequency redistribution is included, a 
good fit to these features can be achieved. Further, to obtain a good fit to the far wings, a small temperature enhancement of the FALF 
model in the photospheric layers is necessary. 
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1. Introduction 

The Second Solar Spectrum is produced by coherent scattering 
processes in the solar atmosphere (Stenflo 1994). Modeling this 
spectrum requires the solution of the polarized radiative transfer 
(RT) equation. It is well known that quantum interference be- 
tween the fine structure (/) levels is responsible for the formation 
of line pairs suchasNaiDi andD2,CanHandK, etc. The linear 
polarization in the Ca n H and K lines was observed and modeled 
by Stenflo (1980), who developed a theoretical framework which 
for the first time demonstrated the profound role that quantu m in- 
terference between different fine structure components within a 
multiplet can have on the observed polarization profiles. An ex- 
tended version of this theoretical framework, applicable to any 
multiplet, was presented in Stenflo (1997). It was used for RT 
modeling of /-state interference in the Nai Di and D2 lines 
by Fluri et al. (2003). The /-state interference theory used in 
the mentioned papers assumed frequency coherent scattering. 
Recently Smitha et al. (2011a, hereafter PI) have extended the 
theory of Stenflo (1997) to include partial frequency redistri- 
bution (PRD) with /-state interference. It is restricted to the 
case of a two-term atom and uses the assumption that the lower 
term is unpolarized. Landi Degl'Innocenti et al. (1997) have pro- 
posed an alternative approach to the same problem based on the 
concept of metalevels, which is intended for dealing with PRD 
in multi-term atoms. A multi-term formulation of /-state inter- 
ference that accounts for the polarization in all the terms, but 
which is restricted to the case of complete frequency redistribu- 
tion (CRD), has been given in Landi Degl'Innocenti & Landolfi 
(2004, hereafter LL04). In the mentioned theoretical formula- 
tions collisional frequency redistribution was not included. 



In Smitha et al. (201 lb, hereafter P2) the redistribution ma- 
trix for the /-state interference derived in PI was incorporated 
into the RT equation, which was solved for simple isother- 
mal model atmospheres. Several theoretical aspects of radia- 
tive transfer in a hypothetical doublet line system are studied in 
that paper. The purpose of the present paper is to perform one- 
dimensional RT modeling of the polarimetric observations of a 
multiplet where /-state interference is relevant. For this we have 
selected the Cn triplet at 5204.50 A (line-1: J b = 1 -» J a = 2), 
5206.04 A (line-2: J b = 2 -> J a = 2), and 5208.42 A (line-3: 
Jb — 3 — > J a — 2). Hyperfine splitting can be neglected because 
the most abundant (90%) isotope of Cn has zero nuclear spin. 

Kleint et al. (2010a,b) have used the Cn triplet for a synop- 
tic program to explore solar cycle variations of the microturbu- 
lent field strength. Recently, Belluzzi & Trujillo Bueno (2011) 
applied the density matrix theory described in LL04 (which is 
based on the CRD approximation) in order to perform a basic 
investigation on the impact of /-state interference in several im- 
portant multiplets in the solar spectrum including also the Cn 
triplet. In this work they have neglected RT and PRD effects. 
However, they have included the effects of lower term polariza- 
tion and the dichroism. They identify and explain qualitatively 
the observational signatures produced by /-state interference in 
the Cr 1 triplet (i.e., the cross-over of Q/I about the continuum 
polarization level occurring between the lines, and the Q/I fea- 
ture around the line-1 core). 

In Section [2] we briefly present the basic equations required 
for realistic RT modeling of lines in the two-term atom picture. 
In Section[3]we present the polarimetric observations of the Cn 
triplet. Section|4]is devoted to a description of realistic modeling 
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of the observations. In Section [5] we present the main results. 
Concluding remarks are given in Section [6] 



2. Polarized line transfer equation for a two-term 
atom 

In a non-magnetic medium, the polarization of the radiation field 
is represented by the Stokes vector (/, Q) T , where positive Q is 
defined to represent linear polarization that is oriented parallel to 
the solar limb. In a medium that is axisymmetric around the ver- 
tical direction, it is advantageous to use a formulation in terms of 
the reduced Stokes vector I instead of the traditional Stokes vec- 
tor (/, Q) T . The transformations between the two can be found in 
Appendix B of Frisch (2007). The relevant line transfer equation 
for the 2-component reduced Stokes vector is 



dI(A,n,z) 



= -k m U,z)[I(A,n,z)-S(A,z)], 



(1) 



in standard notation (see Anusha et al. 201 1). z is the geometric 
height in the atmosphere. See P2 for details of Eq. (Q3 and re- 
lated quantities. The total opacity k lot (A, z) - t]m(A, z)+<t c (A, z) + 
k c (A, z), where cr c and k c are the continuum scattering and con- 
tinuum absorption coefficients, respectively. The line absorption 
coefficient for the entire multiplet is 



T] M (A,Z) = k M (z)4>M(A, Z) =Y J k K j b j a) (l)(Aj h j a ,z), 



(2) 



J ah 



where ki(j b j a ) is the wavelength averaged absorption coefficient 
for the J a — > Jb transition with the corresponding profile func- 
tion denoted by 4>(Aj b j a ,z)- Jb and J a are the total angular mo- 
mentum quantum numbers of the fine-structure levels for the 
upper and lower terms (see Fig. [T}. ki^iz) is the wavelength av- 
eraged absorption coefficient for the entire multiplet. For our 
case of a two-term atom, we need to define a combined profile 
function that determines the shape of the absorption coefficient 
across the whole multiplet. It can be shown that for the Cr i triplet 
line system 0*f(/l, z) is given by (see Eqs. (7) and (8) in Section 2 
ofP2) 



4>mU, z) 



3</>(A l2 ,z) + 5<p(A 22 ,z) + 70(^3 2 , z) 
15 



(3) 



The reduced source vector is defined as 

k M {z)tpM{^ z)Si(A, z) + cr c (A, z)S c (A, z) 



S(A,z) 



ktot(A, z) 
k c (A,z)S th (A,z) 
ktot(Az) 



(4) 



for a two-term atom with an unpolarized lower term. Here <S t h = 
(fix, 0) r , where B^ is the Planck function. The continuum scat- 
tering source vector is 



S c {A,z) 



\L 



$(n')I(A,ix',z)dn', 



(5) 



Since continuum polarization can be seen as representing scat- 
tering in the distant wings of spectral lines (in particular from 
the Lyman series lines, cf. Stenflo 2005), we are justified to use 
the assumption of frequency coherent scattering for the contin- 
uum. The matrix T* is the Rayleigh scattering phase matrix in the 
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Fig. 1. The term diagram showing transitions in the Cr I triplet. 
The diagram is not drawn to scale. 



reduced basis (see Frisch 2007). The line source vector 



f + °° 1 
Si(A,z) = eS ± (*,z)+ - 

J -co ^ 

^ +1 R( A,A',z) 
iU,z) 



x; 



^(n')I(A',iu',z)dfi'dA'. (6) 



The thermalization parameter e = F//(Fr + T/) where Fr and T/ 
are the radiative and inelastic collisional de-excitation rates, re- 
spectively. fi(A, A',z) is a (2 x 2) diagonal matrix with elements 
fi - diag CR ,^), where K K are the redistribution functions 
which include the effects of /-state interference between differ- 
ent line components in a multiplet. K K represents a linear com- 
bination of redistribution functions of type-II and type-Ill. In the 
reduced Stokes vector basis, the angular phase matrix and the 
frequency redistribution functions are decoupled. The phase ma- 
trix part is built into the transfer equation through the T* matrix. 
The theory of redistribution matrices for the /-state interference 
in two-term atoms for the collisionless case is presented in PI. 
This frequency redistribution part that includes /-state interfer- 
ence and the collisional redistribution is given by 



K K (x, x') 



3(2L b + 1) 



25 + 



r Z ^ J, - J ' 



*{rj h r%) JaJf (2J a + l)(2J f + l)(2J b + l)(2J b , + 1) 
L a Lb 1 ) ] L a Lb 1)1 L a L b 1 



J h J f S j \ J b J a S j \ J b ' Jf S 
L a L b 1 1 f 1 1 K 1 f 1 1 K 



Jb' J a S J 1 Jb' Jb Ja f\ Jb' Jb Jf 



(7) 



The ensemble averaged coherency matrix elements (see e.g., 
Bommier & Stenflo 1999) in the above equation are given by 

(r Jb r'j b> )j a j f =AcosPj b ,j b 
x[cospj b , Jb (h\ Jv ) Vf - smB hl j b (fl Jbl )j a j s } 

+B (K) cosBj l/Jb cosa^ Jb 

x{ cos (B Jl/]b + o«J [%(hZ M ) ~ »(/2U/,) 



sin 



(fci^ + affij 



^\ h J h J„J h ,J,\ + *Afj h J a J,.-J. 



(8) 
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The (h l j j ,)j a j f and (fj r ,)j a j f are the auxiliary functions de- 
fined in Eqs. (14) and (15) of PI but are used here for the non- 
magnetic case and with angle-averaged redistribution functions 
of type-II. The auxiliary functions of type-II derived in PI rep- 
resents generalizations of the corresponding quantities appear- 
ing in Sampoorna (2011, see Eqs. (22) and (23)) using a semi- 
classical approach. The important difference between the two in 
the presence of the magnetic field is that in case of /-state inter- 
ference these auxiliary functions depend on both / and m quan- 
tum numbers, unlike in case of ra-state interference where they 
depend only on m. In the particular case of non-magnetic /-state 
interference theory, these quantities depend only on / quantum 
numbers, whereas the corresponding quantities in the m-state in- 
terference theory simply reduce to the well known type-II redis- 
tribution functions of Hummer (1962). Therefore for the nota- 
tional brevity even in the /-state interference case we refer to 
these auxiliary functions as Rn hereafter (in the standard nota- 
tion of Hummer 1962). 

The auxiliary quantities for type-Ill redistribution (see the 
last two lines of Eq. ([8]l) in the non-magnetic case are also 
derived in a way similar to the two-level atom case given in 
Eqs. (27)-(36) of Sampoorna (2011). But in doing so the fol- 
lowing assumptions are made 

1 . Infinitely sharp lower term. 

2. Unpolarized lower term. 

3. Weak radiation field limit (i.e., stimulated emission is ne- 
glected in comparison with spontaneous emission). 

4. Hyperfine structure is neglected. 

5. The (inelastic) collisions that transfer polarization between 
fine structure levels of a given term (upper or lower) are ne- 
glected. 

6. The depolarizing elastic collisions that couple m-states be- 
longing to a given fine structure level Jt, are considered and 
taken to be same for a given term. 

7. We restrict our attention to the linear Zeeman regime 
(Zeeman splitting much smaller than the fine structure split- 
ting). This is not applicable in the Paschen-Back regime. 

The assumption of an unpolarized lower term is made for 
the sake of mathematical simplicity. The inelastic collisions 
that transfer polarization between the fine structure levels are 
neglected. This is justified because the colliding particles are 
isotropically distributed around the radiating atom. This situa- 
tion is similar to the case of an atom immersed in an isotropic 
radiation field producing no atomic polarization (scattering po- 
larization). Neglecting such inelastic collisions is particularly 
valid in the linear Zeeman regime in which we are interested. 
However these inelastic collisions do cause population transfer 
between the fine structure levels, and are properly accounted for 
in the calculations of line opacities (see Section|4]i. 

The angle fij b ,j b is defined in Eq. (10) of PI. The angle Oj ', 
is defined as 



the atomic states. A and B (K) are the branching ratios which are 
given by 



tana 



Jh> Jb 



Uj b ,J„ 



(9) 



Here ha>j b ,j b is the energy difference between the Jy and //, 
states in the absence of a magnetic field. D {K) is the 2K multipole 
depolarizing elastic collisional destruction rate. In general Z) (A:) 
depend on the / quantum numbers of the fine structure states. As 
a simplifying assumption, we take these rates to be the same for 
all the fine structure states of the upper term, and use the classi- 
cal value D (K) = T £ /2 given by Stenflo (1994) where T E is the 
elastic collision rate which is responsible for the broadening of 



r R + r, + r. 



B (K) 



D (K) 



r s + r 7 + £>(*> r« + r / + r £ ' 



(10) 



These branching ratios are the ones derived for a two-level atom 
model by Bommier (1997). In view of the simplifying assump- 
tions stated above, we continue to use the same branching ratios 
for the two-term atom case also. 

The computation of angle-averaged type-Ill redistribution 
functions (that appear in Eq. (O) is very expensive, especially 
for the case of a realistic model atmosphere. Therefore we pre- 
fer to use the approximation of CRD in place of type-Ill redistri- 
bution functions (see Mihalas 1978). We have verified by direct 
numerical computations that this replacement is valid and gives 
results which are almost identical to the explicit use of type-Ill 
redistribution functions. The necessary settings of the branching 
ratios in order to go to the limit of pure CRD are A = and 
B« = (1 - e ) (see also Anusha et al. 201 1). 

3. Observational details 

The Q/I spectra of Cri triplet were observed by Gandorfer 
(2000). In this paper, we present new observations of this triplet 
obtained with the ZIMPOL-3 polarimeter (Ramelli et al. 2010) 
at IRSOL in Switzerland. Fig.|2]shows the observations recorded 
on September 6, 201 1, at the heliographic north pole with the slit 
placed parallel to the limb at /v = 0.15. The polarization modula- 
tion was done with a piezo-elastic modulator (PEM). The spec- 
trograph slit was 60//m wide corresponding to a spatial extent of 
0.5" on the solar image. The CCD covered 190" along the slit. 
The effective pixel resolution in the spatial direction is 4 actual 
pixels wide, due to the grouping of each four pixel rows cov- 
ered by a cylindrical microlense, to allow simultaneous record- 
ing of all four Stokes parameters in the ZIMPOL demodulation 
scheme. The resulting CCD images have 140 such effective pixel 
resolution elements in the spatial direction, each element corre- 
sponding to 1.38", and 1240 pixels in the wavelength direction, 
with one pixel corresponding to 7.84 mA. In Fig.|2]only the cen- 
tral part of the spectral window, corresponding to 1050 pixels 
spanning 8.23 A, is displayed. With the PEM it was possible to 
measure simultaneously one linear and the circular polarization 
component. Measurements of the linear polarization component 
Q/I were alternated with measurements of the U/I component 
through mechanical rotation of the analyzer by 45°. In total we 
accumulated for both the components 2000 exposures of 1 sec- 
ond each. 



4. Modeling of Cri triplet 

To model the Cri triplet the polarized spectrum is calculated 
by a two-stage process described in Holzreuter et al. (2005, 
see also Anusha et al. 2010, 2011). In the first-stage a multi- 
level PRD-capable MALI (Multi-level Approximate Lambda 
Iteration) code of Uitenbroek (2001, referred to as the RH-code) 
solves the statistical equilibrium equation and the unpolarized 
RT equation self-consistently and iteratively. The RH-code is 
used to compute the unpolarized intensities, opacities and the 
collision rates. The angle-averaged redistribution functions of 
Hummer (1962) are used in the RH-code to represent PRD in 
line scattering. In the second stage the opacities and the colli- 
sion rates are kept fixed, and the polarized intensity vector I is 
computed perturbatively by solving the polarized RT equation 
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Fig. 2. CCD image showing the (I,Q/I,U/I,V/I) of the Cri 

triplet. The recording was made on September 6, 201 1, near the 
heliographic north pole at a limb distance defined by // = 0.15. 
The grey scale cuts span a range (from black to white) of 0. 1 % 
in Q/I and U/I, while for V/I the cuts are -0.2 % (black) and 
+0.2 % (white). 

with the redistribution matrices defined in Section |2] which are 
derived for a two-term atom with an unpolarized lower term. 

4.1. Model atom and model atmosphere 

The Cr i atom model is constructed for 14 levels (13 levels of Cr i 
and the ground state of Cm), 11 line transitions, and 13 contin- 
uum transitions. The line components of the 5 S - 5 P triplet of 
Cri are considered under PRD. The values of the various phys- 
ical quantities required to build the atom model are taken from 
the NIST atomic data basqj and the Kurucz data basq^. The data 
for the blend lines are also obtained from the Kurucz data base. 
The photo-ionization cross sections are taken from Bergemann 
& Cescutti (2010). The explicit dependence of these cross sec- 
tions on wavelength is computed under the hydrogenic approxi- 
mation. 

Fig. [3^ shows the temperature structure in some of the stan- 
dard model atmospheres of the Sun - namely FALA, FALC, 
FALF (Fontenla et al. 1993), and FALX (Avrett 1995), which 
we have used in our attempts to fit the observed (/, Q/I) spec- 
tra. Models A, C and F of Fontenla et al. (1993) represent re- 
spectively the supergranular cell center, the average quite Sun, 
and the bright network region in the solar atmosphere. FALX 
represents the coolest model with a chromospheric temperature 
minimum located around 1000 km above the photosphere. Our 
attempts to fit the observed (/, Q/I) spectra using these standard 
models will be discussed in Sect. 5.2. 

We show that a reasonable fit could be obtained only with a 
temperature structure modified with a small enhancement in the 
original temperature structure of the FALF model, in the height 
range of 100 km below the photosphere, extending up to 700 km 
above the photosphere (denoted by FALF). Such an enhance- 



1 www.nist.gov/pml/data/asd.cfm 

2 kurucz.harvard.edu/linelists.html 



ment does not affect the center to limb variation (CLV) of the 
continuum intensity as shown in Fig.[3|3. While such a modifica- 
tion of the temperature structure produces insignificant changes 
in the intensity spectra, Q/I turns out to be quite sensitive to 
these changes in the temperature gradient. 

In order to explore the effect of temperature enhancement in 
a given model atmosphere used for computing line and contin- 
uous spectra, we have performed a simple test (similar to the 
Fig. 3 of Asplund et al. 1999, see our Fig. [3};). It is expected 
that ID model atmospheres (like FALF in our case; or all the 
FAL class of models in general) fit the observed CLV data of 
continuum intensity to a good accuracy. To verify this, we have 
plotted the limb darkening function for the whole range of wave- 
lengths, for different ft values. The theoretical limb darkening 
function fits the corresponding observed data better for fi — > 1 . 
The fit is approximate in the limb positions (say /i = 0.1). In 
Fig. [3J; we also show the theoretical curves computed for the 
FALF model (dotted lines). As can be seen, the limb darkening 
function of FALF does not greatly differ from that of the original 
FALF model atmosphere (a maximum relative difference of 15% 
in the extreme limb). Therefore it justifies a slight modification 
of temperature structure in a given model atmosphere to achieve 
a theoretical fit to the Q/I observations. 



5. Results and discussion 

5.1. Comparison between PRD and CRD 

Fig.|4]shows the comparison between the Q/I profiles computed 
using only CRD (to represent frequency non-coherent redistri- 
bution: solid line), only Rn (dotted line), and a combination of 
Rn and CRD (dashed line) (see Section [2] for the definition of 
CRD). As seen from the figure, the CRD profiles do not pro- 
duce the wing peaks on either side of the line center, which are 
clearly seen in the PRD profiles. Also, the /-state interference 
signatures are more prominent in PRD than in CRD. A good fit 
to the observed Q/I can only be achieved through the use of PRD 
(see Section l5T2l . We have verified that it is impossible to fit the 
observed near wing peaks with CRD alone. The pure Rn mecha- 
nism represents frequency coherent scattering in the line wings, 
the use of which alone also fails to achieve a good fit (since it 
produces too large values of Q/I throughout the wings). We find 
that a proper combination of Rn and CRD is essential to obtain 
a good fit to the observations. This can be seen more clearly in 
Fig. |7J where we present a comparison with the observed Q/I 
profile. It is well known that only such a combination can cor- 
rectly take into account the collisional frequency redistribution. 

Therefore the (/?n,CRD) combination is adopted for the com- 
putations in the present paper. 

The effect of elastic collisions is to cause significant depo- 
larization in the line wings. This can be seen from the dashed 
line in Fig.|4^, which shows that due to elastic collisions the line 
wing amplitudes of Q/I are greatly suppressed with respect to 
the corresponding pure Ru case (dotted line). The fact that the 
Cr i line components are formed in the upper photosphere and 
lower chromosphere makes it necessary to include elastic colli- 
sions in our theoretical modeling. The issue of elastic collisions 
is discussed in some detail in Section[2] 

Fig. [4J3 shows the effect of spectral smearing that needs to 
be applied to the theoretical profiles in order to compare them 
with the observed profiles, which are broadened by the partic- 
ular spectral resolution that was used in the observations. In 
the absence of smearing the Q/I in the line core computed 
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Fig. 3. Panel (a) shows the temperature structure of several stan- 
dard model atmospheres. FALF represents a model with an en- 
hanced temperature structure of the original FALF model. In 
panel (b) we show the center to limb variation of the 'limb dark- 
ening function', I c (p)/I c {n = 1), where I c (fj) is the continuum 
intensity near the Cri triplet. Panel (c) shows the CLV of the 
continuum intensity for all wavelengths covering the violet to 
the IR regions of the spectrum. The observed data are taken from 
Neckel& Labs (1994). 



with Rn and CRD differ significantly (see Fig. |4^). These dif- 
ferences decrease drastically after application of spectral smear- 
ing. Although in isothermal atmospheric models Q/I computed 
with pure Ru and with CRD are very similar in the line center re- 
gion, the same cannot be expected in computations with realistic 
atmospheres. Indeed Q/I computed with PRD shows a double 
peak structure in the line core region with a dip at line center 
(see Holzreuter et al. 2005 for details). The smearing wipes out 
the double-peaked core structure that we see in Fig. 2^. 

5.2. Comparison with observations 

In this section we compare the theoretical Stokes profiles com- 
puted using several standard atmospheric models of the Sun, 
with the observations. Figs.[5]and[6]show the I/I c and Q/I spec- 
tra. From Fig. [5] we can see that the I/I c is not sensitive to the 
choice of the model atmospheres, whereas Q/I is very sensitive. 
The reason for this sensitivity is the angular anisotropy of the 
radiation field, which is different for different atmospheres. 

From Fig. [3]it is clear that the temperature structure of these 
models are considerably different from each other in the line for- 
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Fig. 4. (Q/I) of the Cr i triplet computed with CRD (solid line), R n 
(dotted line), and a combination of CRD and R n (dashed line) for 
yu = 0.15. The positive Q represents linear polarization parallel to the 
solar limb. The thin solid line in Q/I at the 0.025 % level represents the 
continuum polarization. The microturbulent magnetic field B lul b = 0. 
No spectral smearing is applied to the profiles in the panel (a). In the 
line core the Q/I profiles computed with pure Rn and with a combina- 
tion of Rn and CRD nearly coincide (see the insets in the panel (a)). The 
panel (b) shows a comparison between the Q/I profiles computed with 
CRD and with Rn when we also apply a spectral smearing of 80 mA. 
The line types have the same meanings as in panel (a). 
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Fig. 5. Intensity spectra for a choice of model atmospheres. The 
dotted line represents observations and the solid line represents 
the theoretical profiles. The line of sight is represented by fi - 
0.15. 
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Fig. 6. The polarized (Q/I) spectra computed at pi = 0.15 and 
for a choice of model atmospheres. The line types are the same 
as in Fig. The dashed line represents the level of continuum 
polarization. The observations are taken at fi — 0.15. The B tul -b 
values used for the theoretical profiles are given in Table [TJ 



mation region. We find that a modification of the temperature 
structure at certain range in the atmosphere does not significantly 
affect the emergent I/I c profile. However Q/I is quite sensitive to 
such 'modifications' in the temperature structure in the line for- 
mation region. The theoretical profiles (solid lines) in Fig.[6]are 
computed taking into account the effects of microturbulent mag- 
netic fields (Bturb) with an isotropic angular distribution (Stenflo 
1994). Further, the spectral smearing (see Anusha et al. 2010) 
is done using a Gaussian function with FWHM of 80 mA. The 
use of B tur b is essential to obtain correct line center amplitudes 
of Q/I. The values of B tul -b for the three components of the Cr i 
triplet are different. They are chosen to fit the observed line cen- 
ter amplitudes of Q/I using the FALF model. In this way, the 
microturbulent magnetic field strength is only used as a free pa- 
rameter to improve the fit with the observations. We have made 
no attempt to achieve a good fit to the line center amplitudes of 
the (Q/I) profile computed using FALA, FALC and FALX mod- 
els. FALF provides a better fit to the observed Q/I profiles at the 
cores of the three lines and the interference regions in between 
them. However the far wings still remain poorly fitted even by 
the FALF model. To achieve a good fit to the far wing region of 
all the three lines, we found it necessary to enhance (see Fig. [3} 
the temperature in the layers where the far wings are formed. 

Fig. [7] shows a comparison between the profiles of the Cri 
triplet computed with the /-state interference theory (solid line) 
and the observed data (dotted line). This solid line is the same 
as the dashed line in Fig. [4^, except that it now also includes 
the contributions from B tul -b and a spectral smearing of 80 mA 
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Fig. 7. Comparison between the limb {pi = 0.15) observations (dotted 
line, representing the Q/I spectrum of Fig. [2] averaged along the slit) 
and the theoretical profile (solid line). The centers of the 3 lines, the 
PRD peak positions of line-2 and line-3, and the cross-over wavelength 
positions between the lines are marked with vertical lines. The dashed 
line in Q/I at 0.025 % represents the continuum polarization level. The 
solid line is the same as the dashed line of Fig.|4ji, except that we have 
now introduced Hanle depolarization due to microturbulent magnetic 
fields (see Table [TJ, and added spectral smearing of 80 mA to simulate 
the observations. The smearing wipes out the double-peaked core struc- 
ture that we see in Fig. [4k 



to simulate the observations. The best fit values of B tul -b for the 
3 components of the Cr i triplet are given in Table flj The ap- 
proximate heights of formation given in Table [TJ are the heights 
at which the condition t(Aq)/h ^ 1 is satisfied for /i = 0.15. 
The quantity t(Aq) is the total optical depth at line center for the 
FALF model. 

The observed Q/I spectra of the Cr i triplet have two main 
characteristics, namely (i) the presence of a triple peak struc- 
ture in line-2 and line-3; (ii) the cross over of Q/I about the 
continuum polarization level, occurring between the line com- 
ponents. These two aspects are well reproduced in terms of the 
theoretical framework with the redistribution matrix theory for 
/-state interference developed in PI and P2. The small discrep- 
ancies between the observations and theoretical profiles in Q/I 
can be attributed to the presence of blend lines. The blend lines 
are assumed to be formed under LTE and generally depolarize 
the wings of the main line as well as the continuum polariza- 
tion. To fit the observed / spectra the oscillator strengths of the 
blend lines from the Kurucz data base are used unchanged, with 
the single exception of the Y n line at 5205.75 A, since the value 
from the data base for this line does not reproduce the Stokes / 
spectrum at all. Therefore the oscillator strength for this line is 
changed substantially to get a good Stokes / fit. As soon as the 
Stokes / fit becomes good, the Q/I fit automatically becomes 
good as well around 5205.75 A. Such enhancement of the Yn 
line oscillator strength is also used in computing the theoretical 
profiles shown in Figs. [5] and [6] The discrepancy in the theoret- 
ically computed and observed intensity spectra of other blend 
lines is slightly model atmosphere dependent, particularly for 
the one at 5208.6 A (see Fig.[5]for details). We have not made 



H. N. Smitha et al.: ./-state interference signatures in the Second Solar Spectrum 



Table 1. Microturbulent magnetic field strengths necessary to 
obtain a best fit to the line center value of the observed Q/I using 



FALF. 
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a detailed attempt to simultaneously fit the (7, Q/I) spectra of all 
the blend lines. 

The reasons for the lack of a good fit to the observed Q/I 
at the left-wing peaks of line-2 and line-3 remain unclear and 
need further investigation. The deviations of the model fit from 
the observations are possibly due to unidentified opacity sources. 
These deviations however do not affect the diagnostic potential 
of the Cr i triplet. 



6. Conclusions 

In the present paper we have studied the importance of /-state in- 
terference phenomena with realistic radiative transfer modeling 
of the Second Solar Spectrum. We have selected the Cr i triplet 
for this purpose and made use of the theory for partial frequency 
redistribution with /-state interference developed in PI and P2 
in the absence of lower term polarization. This theory is used 
in combination with realistic atmospheres and a model atom for 
Cr i. Our results demonstrate that it is indeed possible to obtain 
an excellent fit to the observed Q/I profiles without use of lower 
term polarization, but they also clearly show that accounting for 
the PRD mechanism is essential to model the observed scattering 
polarization in sufficient detail. The CRD approach alone cannot 
be used to model the observed spectra. We note that Belluzzi 
and Trujillo Bueno (201 1) have carried out a basic investigation 
of the /-state interference phenomenon on different multiplets, 
neglecting RT and PRD effects (the theory they apply is based on 
the CRD assumption). Nevertheless, they were able to identify 
and explain qualitatively the observational signatures produced 
by /-state interference in the Cri triplet (i.e., the cross over of 
Q/I about the continuum level occurring between the lines, and 
the Q/I feature around the line-1 core), neglecting and including 
the effects of lower-term polarization and dichroism. 

Our observations were performed in non-magnetic regions, 
but we find that microturbulent magnetic fields with an isotropic 
angular distribution are needed to fit the line center amplitudes 
of the Q/I spectra. 

The near wing PRD peaks and the characteristic cross-overs 
in Q/I that are typical of /-state interference are well mod- 
eled only through a weighted combination of partially coherent 
(through Rn) and completely non-coherent (through CRD) scat- 
tering processes. The weighting factors (branching ratios) are the 
ones used to represent the collisional frequency redistribution in 
line scattering, and they are properly accounted for in our RT 
calculations. We find that elastic collisions indeed play a ma- 
jor role in modeling the wing polarization of the Cr i triplet. A 
hotter model atmosphere (FALF) with a slight additional tem- 
perature enhancement is found to be needed to obtain a good 
fit to the observed data, in particular for Q/I. This emphasizes 
that the Q/I spectrum (together with the / spectrum) provides a 
much stronger constraint on the model atmosphere than the in- 
tensity spectrum alone. The Second Solar Spectrum is thus not 



only useful for magnetic field diagnostics, but also for modeling 
the thermodynamic structure of the Sun's atmosphere. 
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